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It is assumed that a warm plasma may be described in terms of a continuum theory 
of fluid dynamics in combination with Maxwell's equations. The motions of the heavy 
ions are neglected but collisions with electrons are accounted for by a constant collision 
frequency. For these conditions, a solution is given for the reflection coefficient when a 
plane wave is incident obliquely onto a homogeneous half space of plasma. A d-c magnetic 
field is superimposed on the plasma in a direction parallel to the interface between the 
plasma and the dielectric. It is shown that previous solutions for a cold anisotropic plasma 
and a warm isotropic plasma are recovered as special cases. 



1. Introduction 

In theoretical studies of wave propagation in 
ionized media, it is usually assumed that the plasma 
is cold. This is particularly true in the vast literature 
devoted to ionospheric radio propagation [e.g., Rat- 
cliffe, 1959]. In many cases this is probably well 
justified because the kinetic temperature of the 
plasma is sufficiently low. Nevertheless, it would 
seem to be desirable to investigate the consequences 
of this cold plasma assumption in some specific 
problems. For example, in studying the propagation 
of VLF radio waves, the ionosphere may often be 
approximated by a sharply bounded homogeneous 
plasma [Wait, 1962]. In this paper an explicit 
expression for the reflection coefficient is derived 
under the condition that the electroacoustic velocity 
is finite. In order to simplify the discussion, the 
superimposed d-c magnetic field is taken to be 
parallel to the interface and transverse to the direc- 
tion of propagation. In the limiting case of a cold 
plasma this problem reduces to one first solved 
explicitly by Barber and Crombie [1959]. 

2. Preliminary Considerations 

The plasma medium is regarded as a one-component 
electron fluid. In other words, the ions are neglected 
in the equation of motion, yet their presence is 
required to neutralize the plasma. However, col- 
lisions between electrons and the neutral particles 
are accounted for by an energy independent collision 
frequency v. It is also assumed that the amplitude 
of the plasma and the electromagnetic oscillations 
are sufficiently small that linearized equations are 
valid [Oster, I960]. It is further assumed that the 
drift velocity of the electrons is zero, so the plasma 



as a whole is considered stationary. The average 
number density of electrons is denoted n 0y which is 
regarded as constant in the plasma region. The 

mean velocity is v, while their pressure deviation is 

/;. The electric and magnetic fields are designated 

-» -> 

by their usual symbols E and 77, respectively. 

A uniform magnetic field B is assumed to be 
impressed throughout the plasma. This constant 
magnetic field is taken to be parallel to the y axis 

of a suitably chosen Cartesian coordinate system 

-> -> -* 

(x, y, z). The field quantities E, II, p, and v vary 
witli time t as exp (ico/) where o> is the angular 
frequency. 

For the situation described above, the linearized 
equation of motion for the electrons is 



mno(v+io))v=n Q e(E+vXyB Q )—vp ) 



(1) 



where y is a unit vector in the y-direction, and e and 
m are the charge and mass of the electron, respec- 
tively. As a consequence of assuming a scalar 
pressure term, Landau damping in the plasma is 
neglected. The linearized equation of continuity, 
combined with the equation of state, leads readily to 



u 2 mnoV-v=—iup, 



(2) 



where u is the velocity of sound in the electron gas. 
Maxwell's equations for the plasma in the absence 
of sources are given by 



and 



VXE= — ifx o)H, 
VXlI=ieouE J m ev, 



(3) 

(4) 
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where mo and e are the permeability and dielectric 
constant of free space, respectively. 

In the problem considered in this paper, there are 
two distinct regions separated by a plane interface at 
z=0. The plasma region occupies the upper half 
space 2>0, while a dielectric of electrical constants 
jLi and e d occupies the lower half space z<C0. Follow- 
ing the assumption of earlier workers [Cohen, 1962; 

Hessel et al., 1962; Fedorchenko, 1962], the boundary 

-> 
condition on the electron velocity v is that its normal 
component be zero at the interface 2=0. The ques- 
tion whether an actual interface between a dielectric 
and plasma is totally effective in this regard is not 
answered here. However, the reader is referred to an 
illuminating discussion of this point by Cohen [1962]. 
Very recently the same boundary condition has been 
used by Yildiz [1963], who treated the oscillations of 
an isotropic compressible plasma sphere. 

To simplify the problem, the obliquely incident 
wave from the lower half space propagates transverse 
to the magnetic field. In other words, the wave 
normal is contained in the (x, z) plane and makes an 
angle 6 with the negative 2-axis. Without losing- 
further generality, it is necessary to consider two 
distinct cases. The first and most interesting situa- 
tion is when the magnetic field of the incident wave is 

-> 
parallel to the interface. In other words, H {i) of the 
incident wave has only a y component H { y i} . There- 
fore, we choose it to be of the form 

H { y i) =H e- UQ2 e-^ x y (5) 

where u = i(e d iJLo) l/2 o) cos and A= (e^o) 1/2 co sin 6. 
For the problem as defined, it is found that the 
necessary equations and the boundary conditions are 
satisfied if the reflected wave also has only a y 
component H y r) . In fact, we may write 



H y r) =R l{ H e +u ' z e-^ x } 



(6) 



where, by definition, R u is the reflection coefficient. 
The total magnetic field H y in the region z<^0 may 
be written 

H y =H [e- ikCz +R {l e ikCz ]e- ikSx , (7) 

where (7= cos 6, S=sin 6, and k= {e d ^) l/2 co is the wave 
number in the dielectric region. The electric field 
components in the dielectric are obtained from 
Maxwell's equations. Thus, 



E z 



1 dH y 

i e d co dx 



and 



E x =- 



1 bH v 



ie d a) dz 



(8) 



(9) 



for the region 2<0. The normal wave impedance 
Z, at the interface, is denned by 



7== /Umit\E £ 
\z-^0jH y ' 



(10) 



On combining (7) and (9), it is seen that 



al) 



where r) = k/(e d u) = (iJL /e d ) 1/2 is the characteristic im- 
pedance of the dielectric. Conversely, the reflection 
coefficient may be expressed conveniently in terms 
of the normal surface impedance via the relation 



Rw 



rjC-Z 

'vC+Z 



(12) 



It is apparent from the above development that the 
problem boils down to finding an expression for the 
surface impedance of the compressible magneto- 
plasma at the interface z=0. Because of the con- 
tinuity of the tangential fields E x and H y , this im- 
pedance is simply Z. 

The second case alluded to above corresponds to 
the situation where the electric field of the incident 
wave is parallel to the interface and has only a y 
component. As it turns out, this is a rather trivial 
case from the electromagnetic viewpoint. Neither 
the d-c magnetic field nor the pressure play any role. 
Thus, the solution for the reflection coefficient may 
be obtained by regarding the half space as a cold 
isotropic plasma. This case will not be discussed 
any further. 

3. Development of Equations for the Plasma 

As specified in the preceding section, the magnetic 
field need have only a y component H y . In what 
follows, the subscript y on this quantity is dropped. 
Furthermore, because of the assumed direction of 
incidence d/dy=0. Consequently, (1) in component 
form may be written 



and 



(v+ico)mn v x =n e(E x —v 2 B ) — dp/dx, (13) 
(v J rico)mn v z =n e(E z +v x B )—dp/dz. (14) 



By simple algebraic manipulation, (13) and (14) 
may be rewritten in the form 



gma 



E x 



eo) c 



mcfa 



E z 



1 



dp 



dp 



gamn dx g 2 mn a dz 



(15) 



eoo c 



mg 2 a 



EL- 



E> 



dp 



1 



dp 



gma g 2 mn a dx gamn dz 



(16) 



where a=l-\-(o) c /g) 2 , where g=v-\-iu is a convenient 
complex frequency parameter and o) c = — eB /m is 

the (angular) gyrofrequency for electrons. 
-> 
Maxwell's curl H equation in component form is 

dU 

—^=i€ a)E x +n ev x , (17) 
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and 



dH 
dx 



=i€ Q coE 2 -\-n ev 2 . 



(18) 



This pair may be used in conjunction with (15) and 
(16) to eliminate E x and E z . Thus, 



ieo) c [~dll ie co dp 



mg 2 coe Ka\__dx n e dz 



o 



and 



iugj VdH i€ cc dpi 
: Ka \__dz n Q e dxj 



micogeQ. 



(19) 



ieoo c [~dH ie a) dp 



+ 



mg 2 <joeoKa. |_ dz n e dx 



a 



\ ioog) VdH ^ie Q o) dpi 
in ioo(je {) Ka \_0.r n {) ( dr. J 



(20) 



where o%= (n e 2 ) / (e m) is the (angular) plasma fre- 
quency, 

a=l + (cc c /g) 2 ,K=IQ-Kl 
with 



K l =(l+-r^-\ find K 2 =4 



COqCOo 



Equations (17) and (18), used in conjunction with 
(15) and (16), may be utilized to eliminate v x and v z . 
From this process, we obtain 



E x =- 



K, dH 



it uK dz e uK dx 



K 2 dH Kt-Kbp iK 2 dp 



n eK dx n eK dz 
(21) 



and 



E z - 



K, dH K 2 dH Kr-Kdp iK 2 dp 



ie GoK dx e coK dz 



n eK dz n eK dx 

(22) 



It is evident from the above that all quantities of 
interest are derivable from the two scalars p and H. 
By inserting v x and v z , as given by (19) and (20), into 
(2), we obtain 






v 2 H-[v 2 -r 2 P ]p=o, 



(23) 



where 



dx 2 



d 2 



f-rF-v and Tl 



iwgaK 



( 



1+7 



<4\ 

icogj 



In a similar manner, by working with (3), (21), and 
(22), we arrive at 



where 



( v2 - r ^+f27^=°> 



(241 



Mo€ a>. 



K x 



Equations (23) and (24) are coupled differential 
equations for H and p. An important special case 
is when i? =0. Then, because K 2 =0, these two 
equations decouple. Thus, 



where 

and 
where 



(v 2 -t!)^=o, 

,2 



y} 



_io>9 ( 
u 2 \ 

(V 2 - 7 *)H= 



iwgj 



(25) 



(26) 



-/i «ow 






These are the governing equations for propagation 
in an isotropic compressible plasma when collisions 
are accounted for [Wait, 1964]. 

Another significant special case is when the 
temperature of the plasma is set equal to zero. This 
corresponds to r /; ^oo whence, from (23), the 
pressure p vanishes. Equation (24) then simplifies 
to 

(y 2 -T 2 e )H=0, (27) 

where T 2 e is defined exactly as in (24). This is the 
wave equation for transverse propagation in a cold 
anisotropic plasma [Wait, 1962]. 

The genera] coupled differential equations given 
by (23) and (24), in the case of zero collisions, are 
the same as the ones derived by Seshadri [1963] who 
treated the radiation from a magnetic line source 
embedded in a homogeneous compressible plasma. 

4. Solution of the Coupled Equations 

To obtain the solution of the coupled equations, it 
is assumed that the waves in the plasma vary accord- 
ing to the factor exp (—uz) exp (—i\x) where \=k 
sin 6, and u is to be determined. Because V 2 may 
be replaced by u 2 —X 2 it follows, from (23) and (24), 
that 



(u 2 -\ 2 -T 2 e )(u 2 -\ 2 -T 2 p )-8(u 2 -\ 2 ) 2 =:0, (28) 



where 



iu c K 2 A cogV 1 
gK x V ^ iugj 



When the gyrofrequency approaches zero, the 
coupling factor 8 vanishes, and the two solutions are 
u 2 = \2J r Y 2 e and X 2 +r|. In accordance with earlier 
work, these may be described as electromagnetic 
and acoustic waves, respectively. 
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Equation (28) is now conveniently written in the 
form 



where 

and 

with 



(u 2 -\ 2 -T 2 e )(u 2 -\ 2 -T 2 P )=0, 



T%=v+{tf—w) 112 , 



(29) 



T 2 p =v-(v 2 



n+T 2 e 

"2(1-5) 



and w 



w) 1/2 , 

Tin 



i-d 



The sign of the radicals is chosen so that 1^ 

approaches T 2 e and T% approaches Tp as 6 tends 
uniformly to zero. The solutions of (29) then lead 
to two kinds of waves which are described by the 
factors 

exp (—u e z) exp (—i\x) and exp (—u p z) exp (—i\x), 

where 

^ = ( X 2+r*) 1/2 and u p =(X 2 +T 2 p ) 1/2 . 

It is appropriate to describe these as quasi-electro- 
magnetic and quasi-acoustic waves, respectively. 
The radicals in the expression for u e and u p must be 
chosen to satisfy radiation conditions as £->°o . In 
the case of a finite collision frequency v, this is 
equivalent to assuring that the real parts of u e and 
u p are positive. 



5. Derivation of the Surface Impedance 

General expressions for the magnetic field H and 
the pressure p within the plasma may be expressed 
in the form 



H=[J e exp (—u e z)+J p exp (—u p z)] exp (— i\x), 



and 



(30) 



p=[g e exp (— u e z)+g p exp (— u P z)] exp (—i\x), 

(31) 

where f e , g e , f p , and g p are coefficients which do not 
depend on the coordinates. By using (24), which is 
valid for all 2>0, it follows that 



(32) 
(33) 



and 






A n e 4 

g e =A e — j e > 
e co 


where 






— A n ° e -f 
9p a p ^ ,Jp j 
e co 




A e ~- 


r 2 — 

A 

r 2 

1 e 


A A 

T ' Kl and A _ r '~ r ^i 
K 2 ana A » * K 2 

1 p 



Thus, the pressure is given in terms oiJ e and/^ by 
P^— [f e A e exp (—u e z) +J p A p exp (—u p z)] exp (—i\x) . 

€()C0 

(34) 

Using (20), (30), and (34) and the boundary 
condition v z =0, we obtain 



/, 



Q a Jg)[u.-A*\-(l+j£)(k-A, 



U e ) 



f. 



(io> e /g)[u P 



-A p \]-(l+S^(\-A P u P ) 



(35) 



This means that the rigidity condition at the interface 
fixes the ratio of the quasi-acoustic wave to the 
quasi-electromagnetic wave. In the special case of 
a cold plasma (i.e., u=0), A^-^oo and, therefore, 
the ratio j v \j e vanishes. Another limiting case is 
when the plasma becomes isotropic (i.e., co c =0). 
Then, again, A p ->oo and j p \j e become vanishingly 
small. However, in this case, a more meaningful 
ratio is g p /fe, which approaches the finite value 
noe\/(e a)U p ) in accordance with the required be- 
havior for an isotropic compressible plasma [Wait, 
1964]. 

Using (21), (30), and (34), the following expression 
is obtained for the tangential electric field at the 
interface: 



E T 






e co 
e cx) |_ iK 



K 2 i\ , Kj-K . .. K 2 . . 
— jFrH ^ — A^X — jy A e %u R 



K ' K ^ e ° ,x K 1VeljWe y e 

K 2 i\ x K x -K . . K 2 . . 1, 

J\ p %h jr A-ptllp \Jp. 

(36) 



K 



K 



According to (10), it is noted that 

E x ] z = =ZH v ] z =o=Z(f e -\-J p ) . 



(37) 



Therefore, on equating (36) and (37), we arrive at 
the rather cumbersome expression for the surface 
impedance. 



<- 



\4) 



K 2 



K 

K x u e 



-(#!-£) A. 



K 2 A e u e 



+ f! K '- 



K x u p 



-(Ki-KDAp- 



X 

K 2 A p Up 



] 



(38) 



where fp/f e is given explicitly by (35). 

As a partial check on the correctness of (38), we 
consider the case of a cold anisotropic plasma (i.e., 
u=0). Then, since fp/j e —>0 and A e -^0, it imme- 
diately follows that Z^Z C where 

Z C =^^{K 2 \-K l u e ) with ^=(\ 2 +H)*, 

€qC0XV 
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which is in agreement with a result derived originally 
by Barber and Crombie [1959]. Another special 
case is a warm isotropic plasma (i.e., co c =0). Then, 

6=0, T e ^T e , r p -*T p , Kz-^0, and K/Ki-+[l + <dil 
(icog)] = e/e where e is the dielectric constant of the 
isotropic plasma, In this limiting case, (f P /fe)->0 
and A p — >oo but in such a way that (A p f p /J e )-^\/Up. 
As a result of these considerations, Z^Z t where 



where 



with 



and 



leco L u„u e gw J 



u e ={\ 2 +yl)i and u p =(\ 2 +yl)$ 



(39) 









7 2 e=— e Mo^ 2 






e/z c<r 



The above result for Z 7 is in agreement with an 
independent derivation for an isotropic warm plasma 
[Wait, 1964]. It is immediately evident that the 
bracketed term in (39) may be replaced by unity 
if u/c<^<^\ , which is the condition for the validity 
of the cold plasma assumption. 

6. Concluding Remarks 

An explicit expression for the reflection coefficient 
for a compressible magnetoplasma has been given in 
this paper. It is shown that it reduces to the well- 



known result for a cold magnetoplasma when the 
ratio of the acoustic velocity u is small compared 
with the velocity c of electromagnetic waves in free 
space. The general expression also reduces to a 
previously derived result for reflection from an 
isotropic compressible half space of plasma. 
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